THE HAUSDORFF DIMENSION OF THE 
BOUNDARY OF THE LEVY DRAGON 



P. DUVALL, J. KEESLING 

Abstract. A theoretical approach to computing the Hausdorff dimension of the topological 
boundary of attractors of iterated function systems is developed. The curve known as the Levy 
Dragon is then studied in detail and the Hausdorff dimension of its boundary is computed using 
the theory developed. The actual computation is a complicated procedure. It involves a great 
deal of combinatorial topology as well as determining the structure and certain eigenvalues of 
a 752 x 752 matrix. Perron-Frobenius theory plays an important role in analyzing this matrix. 



§1. Introduction 

In [LI], P. Levy studied the complex curve which has come to be known as the Levy 
Dragon. He showed that the Dragon is self-similar, and argued that the plane can be tiled 
by copies of it. We give a self-contained account of the Dragon that determines the properties 
of this set that we need in our computations. However, the paper by Levy is a tour de force 
on this subject. It is amazing that he was able to determine so many properties of this and 
other curves in [LI] without the use of modern computers. The reader is urged to see [L2] 
for an excellent annotated translation of [LI]. 
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Figure 1. The Levy Dragon 



Although Levy found so many properties of the Levy Dragon, he did not compute the 
Hausdorff dimension of its boundary. In [L2, page 236] G. Edgar asks whether the dimension 
of the boundary of the Dragon is greater than one. In this paper we answer that question 
in the affirmative and in fact determine a precise value for this dimension. In the process 
of doing this we develop a general theory for computing the Hausdorff dimension of the 
boundary of self-similar fractal tiles. This depends heavily on the paper by K. Falconer [Fl] 
who gave an (abstract) method for determining the Hausdorff dimension of sub-self-similar 
sets. Our contribution is to show that in the case of boundaries, it is possible to derive 
a method which can be carried to a computation. We then make a detailed study of the 
Levy Dragon and show how to apply our theory to compute the dimension of its boundary. 
In subsequent work the authors together with A. Vince have found a general method for 
finding the dimension of the boundary of self-similar digit tiles [DKV]. The size of the 
matrices involved in computing the dimensions of the boundaries of these tiles is typically 
less than 10 x 10. Unfortunately the Levy Dragon is not a self-similar digit tile and so these 
methods do not apply directly to this example. Fortunately, motivated by some ideas in this 
paper R. Strichartz and Y. Wang [SW] found a method of computing the dimension of the 
boundary of the Levy Dragon which has considerably reduced the computations required 
in this case as well. In their method one only needs to determine a certain eigenvalue of an 
11 x 11 matrix. Strichartz and Wang made use of the graph directed constructions of D. 
Mauldin and S. Williams [MW] rather than the sub-self-similar sets of Falconer [Fl] used 
in this paper and in [DKV]. 

First, we need to establish some terminology and notation. We assume that the reader is 
familiar with the basic theory of self-similar fractals. The books by Edgar [El] and Falconer 
[F2] are excellent references. The book [E2] contains a number of papers which give a nice 
historical perspective on the subject. 

If X is a complete metric space and F = {fi, f%, • ■ ■ , f n } is a collection of contraction 
mappings of X to itself, then F is said to be a (hyperbolic) Iterated Function System or 
IFS. It is well known [H] that for such an F there is a unique compact set K such that 
K = (JlLi fi(K)- K i s caue d the attractor or invariant set for F. In fact, if C(X) denotes 
the collection of non-empty compact subsets of X with the Hausdorff metric, then the map 
F : C(X) — > C(X) given by F(A) = \J™ =1 fi(A) is a contraction mapping. Thus K is 
the unique fixed point of F guaranteed by Banach's contraction mapping principle, and 
K = lim F k (A) for any A G C(X). In this paper, we will only be concerned with the case 

k^oo 

in which the fi are contracting similitudes, maps satisfying d(fi(x), fi(y)) = Cid(x,y), where 
d is the metric on X and < q < 1. 

Our main interest in this paper is the Levy Dragon. It is the attractor for the IFS 
F = {fi, f 2 }, where fi and fi are similitudes of Euclidean space R 2 given by 



fi(x,y) 



h(x,y) 
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d(I, J) 



However, much of the theory developed in the early sections holds in a much broader setting. 
So, it is appropriate to have fairly general notational conventions. The Hausdorff dimension 
of a set E is denoted by dim# E, and the lower and upper box-counting dimensions of E are 
denoted by dim p E and dim^E 1 . Definitions of of these dimension functions can be found 
in [Fl], [F2] and [E]. 

Let O denote the collection of sequences i" = {zi, z 2 , ■ • ■ } with ik G {1, 2, • ■ • , n}. will 
denote the set of sequences of length k with entries from {1, 2, • • • , n}. O has the natural 
metric given by 

°ii c i 2 ---Cim! if ^ = 3k for k < m and i m+1 ^ j m+1 . 

For I e fl, we let Ik denote the finite sequence {h,^,- -- Ak}, and will often use the 
shorthand fi k for the composition fi 1 fi 2 ■ • ■ fi k and c\ k for the product Ci 1 Ci 2 ■ ■ ■ Ci k . We 
will also need the maps a : O — > O and g : O — > iiT given by ct(h, Z2, - • •) = {^2, 23, • • • } and 
g(I) = lim fi 1 fi 2 ■ ■ ■ fi k (x). It is not difficult to see that g(I) is independent of the choice 

k — >oo 

of x G X. The map a is the well-known shift map on the sequence space. 

§2. Sub-Self-Similar Sets 

In [Fl], Falconer introduced the notion of sub-self-similar (s.s.s.) sets, and gave an 
(abstract) way to determine the dimension of such sets. In this section, we give a brief 
summary, in our notation, of the notions from [Fl] that we will need for our calculations. 
Given a set F = {/1, fa, • • ■ , f n } of contracting similitudes on R z , the closed set E is said 
to be sub-self-similar for F if E C \J™ = i fi(E). It follows that E C K, where K is the 
attractor for F. Note in particular that the topological boundary dK of K is s.s.s., since 
the interior int K of K is mapped into itself by each of the open mappings /j. The first 
crucial observation about s.s.s. sets is 

2.1 Proposition. Let E be a closed set. Then E is s.s.s. for F if and only if E = g(A) 
for some compact set A C O such that a (A) C A. Furthermore, if E is s.s.s., such a set A 
is given by A = {I G Q\g(a k (I)) G E for all k > 0}. 

Proof. This is Proposition 2.1 (and its proof) in [Fl]. □ 

The set A in Proposition 2.1 is fundamental in the calculation of dim// E. For k > 0, let 
A k be the set of finite sequences obtained by truncating elements of A after k terms. For 
s > define 

i/k 

T ^ = & ( £ ' 

Recall that F satisfies the open set condition [H] if there is a bounded nonempty open set 
U such that 

n 

\Jfi(U)cU,fmd 

i=l 

f i (u)nf j (u) = ib,i£i^j. 
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Our main tool for calculating dimensions is the following theorem of Falconer. 

2.2 Theorem. Let E be a closed set which is s.s.s. with respect to the family {/i, fi, • • • , fn) 
of contracting similitudes which satisfies the open set condition. Let s be the number satis- 
fying t(s) = 1. Then s = dim# E = dim p E = dim^i?. 

Proof. This is Theorem 3.5 of [Fl]. The packing dimension dimp(X) is currently an impor- 
tant dimension function as well. Since dimuiE) < dirrip(E) < dimpiE) is always satisfied, 
we also have that s = dimp(E) as well. □ 

§3. Boundaries of Attractors 

In this section, we specialize to the case of the boundary of the attractor of an IFS satis- 
fying the open set condition, and show that stronger invariance properties of the boundary 
give a more tractable version of Theorem 2.2 for computing its Hausdorff dimension. Let 
F = {fi, f2, ■ • ■ , f n } be an IFS satisfying the open set condition, let K be its attractor, 
and let A C O be as in the previous section. One obvious difficulty in applying Theorem 
2.2 is that one has to sum over finite sequences determined by an infinite condition. Let 

A' = {I e n\fi n (K) ndK^dl for all n > 0}. 

We will show that A' = A and use that to show that membership in A n can be determined 
from finite data. 

3.1 Proposition. a(A') C A'. 

Proof. Suppose that I C A' but a(I) £ A'. Then for some n, 

fiJ i3 ..-f in (K)cmtK. 
But since each /j maps the interior of K into itself, 

f^U 2 - ■ ■ U n {K) C f^{mt K) dint K, 

a contradiction. □ 

3.2 Proposition. A = A'. 

Proof. Suppose that I G A' . Let k be given, and use p G K to compute g. For any m > 0, 
3.1 implies that fi k fi k+1 ■ • • fi k+m {K) fl dK ^ 0, so 

d(fi k fi k+1 ' ' ' fi k+m (p) ) dK) < Ci k Ci k+1 ■ ■ ■ Ci k+m . 

Therefore lim fi k fi k+1 ■ ■ ■ fi k + m (p) G dK and I G A and A' C A. 

m — >oo 

Suppose that I G A, but f h f i2 ■ ■ ■ f ik (K) n dK = 0. Then 

f h f i2 ■■■f ik (K)cmtK. 

Since 

fix /i 2 ' ' ' fik \K) ~> fh /*2 ' ' ' fik + i \K) ^ ' ' ' fil /*2 ' ' ' fik+ m (-^) 

for all m, it follows that g(I) G int a contradiction. Thus A C A' . □ 
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3.3 Proposition. A k = {I e fl k \f Ik (K) n dK ^ 0}. 

Proof. Let I k be a sequence in O fc such that fi k (K) n <9-ftT 7^ 0. If f Ikl (K) C int X, 
we would have that fi k (K) = f Ikl {f ik {K)) C fj k _ 1 (K) C int if, a contradiction. Thus 
// fc _ 1 (iT) fl if 7^ 0, and, working backwards, we get that //. (if) fl if 7^ for j < fc. Next, we 

need to show that ifc extends to a sequence in A. Since // fc = [J/i fc (/j(-K')), there is some 

j 

ifc+i such that • • • fi k fi k+1 (K) H <9if 7^ 0. Continuing by induction gives the necessary 
extension of I k . □ 

Proposition 3.3 leads to a strategy for computing dim# <9if . Suppose for convenience 
that the q all have the same value c, and suppose that we can determine a = lim lAfcl 1 ^, 

k^oo 

where \S\ denotes the cardinality of the set S. Then 




t(s)= lim [ V c? fe ) =lim (\A k \c k f /k 



k— >oo \ z — ' fe / 



ctc s . 



Thus dim# <9if is the solution to ac s = 1, or 

dim H dK 



ln(a) 



ln(c) ' 

§4. The Levy Dragon 

For the rest of this paper, K will denote the Levy Dragon, which is the attractor of the 
IFS F = {/1, ^2}, where fi and f2 are the similitudes of Euclidean space R 2 given by 

, / x ,x-y x + y 

h(x,y) = (— , — g— )> 

y - x + 1 
2/) = ( g ' 2 ^ 

Levy [LI] studied K extensively, and showed, among other things, that K tiles the plane 
in the sense that the plane can be written as the union of congruent copies of K that meet 
only in their boundaries. It is first important that we establish that K satisfies the open 
set condition so that we can use Theorem 2.2. The following Lemma can be found in [B] or 
[K]. 

4.1 Lemma. Suppose that {fi,---,f n } is an IFS on R/ with K the invariant set. Suppose 
that J27=i c f = 1- Then if K has nonempty interior in R/ , then K satisfies the open sets 
condition with the int K being the open set. 

Proof. Let U = int K. Then since each fi is an open mapping in R l , it must be the case 
that fi( int K) C int K for each 1 < i < I. By Lemma 1.1 of [K] it must be the case 
that H l (f % (K) n fj(K)) = for all i ^ j, where H a (A) denotes the Hausdorff ct-measure of 
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A. Since Lebesgue /-dimensional measure is proportional to TL , this implies that A(/i(int 
K) n (/j(int K)) = 0. Of course, this is only possible if ^(int K) n (/j(int K) = 0. □ 

For our study, we will need some more notation. Let S be the unit square with vertices 
{(0,0), (0,1), (1, 1), (1,0)} and let T C S be the triangle with vertices {(0, 0), (1/2,1/2), (1,0)}. 
We will view K as lim F k (T ). 

k^oo 




Figure 2. T And Its Neighborhood 



To and its neighborhood in Tq are shown in Figure 2. S is the union of the four triangles in 
S with hypotenuse an edge of S and third vertex at (1/2, 1/2). We denote the triangulation 
of R 2 consisting of these triangles and their integral translates by T . If T is a right triangle 
in the plane, we refer to the vertex opposite the hypotenuse as the top vertex of T, and the 
other vertices as left and right so that we encounter the vertices in the order left, top, and 
right as we traverse the boundary in clockwise order. The left(right) edge of T is the edge 
determined by the left(right) and top vertices of T. The triangulation T has a subdivision Ti 
obtained by subdividing each triangle T into the two triangles determined by the left(right) 
edge of T and the midpoint of the hypotenuse of T. Continuing in the obvious way gives a 
sequence Tq y T\ >- • • • >- >- • • • of subdivisions of R 2 into right isosceles triangles. Each 
triangle in has diameter (-^) k - For each T G 7^, let T L be the triangle in 7fc + i whose 

hypotenuse is the left edge of T and which has a vertex outside of T. Define T R similarly. 
One can easily see that F(T ) = T L [JTq 1 . In fact, the following Proposition is immediate. 

4.2 Proposition. For each k > 0, F k (T ) is a union of 2 k triangles in T^. F k+1 (T ) is 
obtained from F k (T ) by replacing each T in F k (T ) by T L and T R . □ 

We would like to use the formula at the end of §3 to compute the Hausdorff dimension 
of dK. The problem is that K and dK are very complicated objects, and we can only see 
finite approximations to them, so it is not clear how to compute A^. What we can hope to 
compute is Bt = {Ik £ ^k\fi k (To) H dF k (T ) ^ 0}. We now attack that problem together 
with the question of relating Ak to Bk- 

For T G T fc , let N(T) be the star of T, that is, N(T) = {T G T k \Tf)T' ^ 0}. Each N(T) 
consists of 15 triangles. For technical reasons to be seen later, we order these triangles as 
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N(T)[1],N(T)[2},--- ,AT(T)[15] as follows: T = N(T)[1}; then N(T)[2] - N(T)[8] are the 
triangles incident with the left vertex of T in clockwise order, so that N(T)[8] has the left 
edge of T for its right edge; N(T) [9] is incident with the top vertex of T and and has the 
left edge of N(T)[8] for its right edge;iV(T) [10] has the right edge of T for its left edge; and 
N(T) [11] — N(T) [15] are the remaining triangles incident with the right vertex of T ordered 
clockwise. This ordering is pictured in Figure 3. The neighborhood type p(T) is the 15-tuple 
P( T ) = {Pi; P2, ■ ■ ■ , Pis}, where 



We say that T E T k is of type p(T) and that T is covered if p(T) is the all ones vector. T 
being covered means that T and all of the triangles in T k that intersect T are contained in 
F h (To). To avoid even further notation, we will use N(T) to denote both the collection of 
triangles as defined above as well as the set which is the union of those triangles (the carrier 



One of the difficulties in dealing with K is that it is not easy to locate points in the interior 
of K by looking at F k (T ). One needs to find objects of positive area which persist from 
iteration to the next, but as soon as a triangle T appears at one level, its interior is discarded 
at the next level and is replaced by two triangles exterior to T. It is not until the iterations 
have become sufficiently intertwined to create covered triangles that the discarded mass is 
filled back in by adjacent triangles. Through computer experimentation we discovered that 
covered triangles do not exist until the 14th iteration, and only 8 of the 2 14 triangles that 
comprise F 14 (To) are covered! However, if a triangle T e T k is covered, then T is the union 
of two triangles in T^+i which are also covered. It follows that int T C F k+m (To) for all 
m > 0, so that int T C int K. We record this and a few other elementary observations 
about neighborhoods in the next Proposition. 





otherwise. 



of N(T)). 




Figure 3. The Neighborhood Ordering 
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Figure 4. iV 5 And F 5 T 

4.3 Proposition. 

(1) Let N = AT (T ). Then F(N ) C N , and therefore K C N . 

(2) IfT eT k and T'cT /or some T' G T fc+ i, i/ien N(T') C N(T). 

(3) If T ETk is covered and f CT /or some T" G T/c+i, £/ien T" is covered. 

(4) If T ETk is covered, then intT C intK . 

(5) is a vertex of any triangle fi k (T ) C F fc (T ), it is a/so a vertex of some triangle 
fi k+1 {T ) C F fc+1 (T ), and therefore v E K. □ 

As we mentioned above, Levy [LI] observed that the Dragon iiT tiles the plane, and 
the tiles occur in the following way. For each triangle T a E T there is an IFS F a whose 
similarities are related to T a in the same way that /i and / 2 are related to T . At each 
stage, the F k (T a ) tile the plane, and the limits {K a } give a tiling by copies of the Dragon. 

Now recall our definitions B k = {I k E n k \f Ik (T ) n dF k (T ) ^ 0} and A k = {I k E 

n k \f Ik (K)ndK^®}. 

4.4 Proposition. B k C A k . 

Proof. Note that B k is the set of I k such that fi k (To) is not covered. Suppose that I k E B k , 
and let T = fi k (T ). Then for some T E N(T), T £ F k (T ) and thus V C F k (T a ) for 
some other tile at the kth level. Let v be a vertex common to T and T' . By 4.3 (5) applied 
to both F and F a , v E K n -KT a , so u G <9-ftT. Since n = fi k (w) for some vertex to of T , 
i; G / 7fe (K), so / 7fc (K) n dK £ 0, and 7 fc G A fc . 

It follows that \B k \ < \A k \. By using B n in two ways, we can compute dim# dK. Let 
P = lim |S fc 

k — >0O 

4.5 Theorem 



l/k 



dim H dK = 



HP) 
HV2)' 



Moreover, < H a (dK) < +oo where H a (A) is the Hausdorff a -dimensional measure 



and a 
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Proof. Since f3 < a, we get 

i^L<dim H 3K 

MV2) - 

from the formula at the end of §3. On the other hand, we can also use f3 to get an upper 
bound for the box counting dimension of dK, which is equal to the Hausdorff dimension by 
Theorem 2.2. Let N k denote the the set of triangles T in the subdivision of iVo induced by 
7fc. For such a T, it is clear that if p(T) is the all zeros vector, then int T n K = and if T 
is covered, int T n dK = 0. Thus, dK C |J{^ e N k\p(T) is not constant}. By expanding 
to neighborhoods, we get dK C \J{N(T)\N(T) contains some fi k (T ), G -Bfc}. There are 
at most 15|Sfc| such N(T), so for each k, is covered by at most 15|Sfc| sets of diameter 
< 3(v / 2)~ A: . Therefore, the box counting dimension of dK is less than or equal to 

Um ln(15|S fc |) _ HP) 



fc-oo ln(3(v / 2)- fc ) ln(v^) 

The referee pointed out the second part of the theorem. It is one of the advantages of the 
approach of Strichartz and Wang [SW] using the graph directed construction of Mauldin 
and Williams [MW] that one can claim that < H a {dK) < +oo. It also follows from [MW] 
and the fact in this paper that the matrix M used in arriving at the asymptotic estimate 
of \Bk\ is irreducible. □ 

§5. Computational issues 

In the light of Theorem 4.5, our goal is to understand the asymptotic behavior of the 
number \Bk\ of triangles in F h (T ) which are not covered. To this end, we need to study 
the neighborhood structures introduced in §4. We will index the binary vectors of length 
15 by integers between and 2 15 — 1 = 32767 by the correspondence (xi,X2, ■ ■ ■ ,^15) 
Yll=i Xi2 l ~ 1 . Thus we could indicate that T is covered by writing p{T) = 32767. Recall 
that each T e Tk is the union of two triangles T\ and T2 in Tk+\i where we let T\ be the 
triangle that contains the left vertex of T, and T<i contains the right vertex. If we know the 
neighborhood type for T, it is simple (but tedious) to write down the neighborhood types 
for Ti and T 2 . We have 

5.1 Proposition. For T eT k , if 

p(T) = (xl, x2, x3, x4, xb, x6, x7, x8, x9, xlO, xll, xl2, x!3, x!4, xlb), 



then 
and 



p{Ti) = (x8, xl, x9, x8, xlO, x9, xl, xlO, xlb, x3, x2, xb, x4, xl, x6), 



p(T2) = (xlO, xl, xl2, xll, xl4, xl3, x2, xlb, x3, x8, xl, x9, x8, xlO, x9). □ 
Proposition 5.1 gives a way to compute \Bk\ in principle. Let V(k) be the 2 15 long vector 



such that V(k)i = number of T G N k such that p(T) = i. Then let M be the 2 15 x 2 15 
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matrix with Mjj = t, if a triangle of type i gives rise to t subtriangles with type j. For 
example, M , = 2, M 327 67,32767 = 2, M 0)J = 0, % > 0,and M 32767J = 0,j < 32767. M can be 
constructed from the information in 5.1. Then M has the property that V(k + 1) = V(k)-M, 
so if J is the column vector with Jj = 1 if i is odd and less than 32767, otherwise, then 
\Bk\ = V(0) ■ M k ■ J and we could study the asymptotic behavior of \Bk\ via the eigenvalues 
of M. Fortunately, we do not have to go to such extremes. A little thought shows that not 
every one of the 2 15 neighborhood types can occur. Our computations determined that, in 
fact, relatively few types actually do occur. 

The only triangle in F°(Tq) is To itself, so the neighborhood types for the fifteen trian- 
gles in N are encoded by 1, 2, 2 2 , • • • , 2 14 . For example, T is triangle 1 in its neighbor- 
hood, so p(Tq) = 1. For a more revealing example, consider the triangle T with vertices 
{(0,0), (0,1), (1/2, 1/2)}. It is number 8 in the neighborhood of T , but T is number 
10 in N(T). Therefore, p(T) = 2( 10 " 1 ) = 512. If we let S k = {p(T)\T E N k }, then, 
So = {2 fc ,0 < k < 14}. By applying the formulas above for p(T\) and p(72), we can 
compute Si, S2, and so on. Again, we determined by computation that S19 = S20, so that 
Sk+iQ = S19 for all k > 0. Let be this common "stable" set of neighborhood types, 
Soo = {si, S2, • • • , 5752}- (We will discuss the order on below.) The cardinality of Soo 
is 752, so we can perform the analysis indicated above with state vector and transition 
matrix of manageable size. For k > 19, let V(k) be the vector of length 752 whose ith 
entry is the number of triangles in Nk which have type S{. Let M be the 752 x 752 matrix 
with Mi j = t, if a triangle of type Si gives rise to t subtriangles with type Sj. Then as 
above, V(k + 1) = V(k) ■ M, and if J is the column vector whose ith entry is 1 if is coded 
by an odd number 7^ 32767, and otherwise, then 

\B k \ = V(k)-M- J, 

so that 

|5i 9+fc | =V(19)-M h - J. 

The types in have a nice structure. The types and 32767 are absorbing in that each 
gives rise to two triangles of the same type at the next level. Let S a = {0, 32767}; let 

S t = {4, 8, 9, 16, 32, 64, 66, 128, 512, 1024, 1026, 2048, 4096, 8192, 8193, 16384}; 

and let S e = \ (S a [j St)- It can be shown by computation that if we order the types in 
Soo so that Si < Sj < Sk H Si & St,Sj & S e ,Sk & S a , then M has the block form 

(P Q R 
M= C L 
\ I 

where C is a 734x734 matrix and lis the 2x2 identity matrix. Write V(k) = (V t (k),V e (k),V a {k)), 
and let 

' (?) 
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where the blocking is consistent with that of M. Then we have 

\B 1Q+k \ = (V t (k),V e (k))- 



P Q\ h [Jt 



C ) \ J, 



' e 



The matrix P can be seen to be a permutation matrix, and C is primitive in the sense 
of Perron- Frobenius theory [BP], i.e., all entries of some power of C are positive. In fact, a 
calculation shows that C 25 3> 0. It follows [BP, Theorem 2.4.1] that C has an eigenvalue A 

equal to its spectral radius, and that lim —rC k = D, where D is a positive matrix. In our 

k^oo A fc 

case, A > 1, so 

1 (P Q\ k _ v _ f±P k ^P k ~ 1 Q\_fO 



Therefore, 



Then we have 



By §4, 



^^^ = A^ K(19) ' D ' Je = 9>0 - 



(3= lim \B rn \ 1 / m = lim Xq 1/m = A. 

m^oo m^oo 

ln(A) 



dim# OK 



Floating point computations using the power method give the estimate A ~ 1.954776399. 
That approximation for A gives dirnj^ dK m 1.934007183. According to [BP, Theorem 
2.2.35], the spectral radius r of an irreducible matrix L satisfies u < r < U ', where u 
and U are the minimum and maximum row sums of L, so it follows that our A satisfies 
u k < A < U k for each k, where and Uk are the minimum and maximum row sums 
of C k . We computed u 30 = 67936360 and U 30 = 727212953 using integer arithmetic, which 
gives the rigorous estimate 1.824190 < A < 1.974189. Thus A is provably greater than 
and the Hausdorff dimension of dK is greater than 1. Thus, the boundary of the Dragon has 
Hausdorff dimension greater than 1, and Edgar's question is answered in the affirmative. 
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